Systematic Microcanonical Analyses of Polymer Adsorption Transitions 
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We investigate the cooperative effects of a single finite chain of monomers near an attractive 
substrate by first constructing a conformational pseudo-phase diagram based on the thermal fluc- 
tuations of energetic and structural quantities. Then, the adsorption transition is analyzed in more 
detail. This is conveniently done by a microcanonical analysis of densities of states obtained by 
extensive multicanonical Monte Carlo simulations. For short chains and strong surface attraction, 
the microcanonical entropy turns out to be a convex function of energy in the transition regime. 
This is a characteristic physical effect and deserves a careful consideration in analyses of cooperative 
macrostate transitions in finite systems. 



I. INTRODUCTION 

The understanding of the adsorption phenomena of 
polymers on surfaces is a prerequisite, e.g., for design- 
ing micro- or nanostructures. Also the fact that various 
polymers are usually found near both sides of cell mem- 
branes and are important for their mechanical stability 
and physiological function [l| has driven studies of poly- 
mers near surfaces and interfaces [2]. In this context 
a deeper understanding of the origin of specific binding 
affinities of proteins regarding the type of the substrate 
and the amino acid sequence is very desirable. In recent 
years, some progress has been made in this field [3, 4], but 
due to the complexity introduced by the huge amount 
of possible sequences and surfaces many problems are 
still open. A qualitative understanding of the coopera- 
tive nature of the adsorption transition for short chains 
can, however, already be gained by studying the behav- 
ior of homopolymers close to a flat substrate [Bl-Q • We 
focus here on the systematic description of the phase di- 
agram in a wide parameter range and cooperative effects 
of chains of finite length. 

First, in this numerical study, the conformational 
pseudo-phase diagram of a coarse-grained non-grafted 
off-lattice polymer will be constructed versus temper- 
ature and surface attraction strength. The competi- 
tion between monomer-monomer and surface-monomer 
attraction gives rise to a variety of different conforma- 
tional phases [H, @. Our computer simulations rely on 
the multicanonical Monte Carlo method [To| that allows 
for the precise determination of the canonical expecta- 
tion values of suitable observables over a wide range of 
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temperatures within a single long simulation run. In ad- 
dition it yields an estimate of the density of states, which 
possesses a convex regime at the adsorption transition in 
the case of short chains and strong surface attraction. 
Albeit known to be a continuous transition in the ther- 
modynamic limit of infinitely long chains [5], the adsorp- 
tion transition of non-grafted finite-length polymers thus 
exhibits a clear signature of a first-order-like transition, 
with coexisting phases of adsorbed and desorbed confor- 
mations. 



II. MODEL AND SIMULATION 

Off-Lattice Homopolymer with Attractive Substrate. 
We employ a coarse-grained off-lattice model for ho- 
mopolymers that has also been generalized for studies 
of heteropolymers [11] and helped to understand protein 
folding channels from a mesoscopic perspective [12]. Ad- 
jacent monomers are connected by rigid bonds of unity 
length, but bond and torsional angles are free to rotate. 
The energy function consists of three terms, 
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where the first two terms give the energy of a polymer in 
bulk (£?buik) that consists of the standard 12-6 Lennard- 
Jones (LJ) potential and a weak bending energy. The 
bending energy provides a penalty for successive bonds 
deviating from a straight arrangement. Here < < tt 
denotes the bending angle between monomers z, and 
i + 2. The distance between the monomers i and j is 
and Zi is the distance of the zth monomer to the substrate. 
The third term is the attractive surface potential i? S urf, 
obtained by integrating over the continuous half-space 



z < 0, where every space element interacts with a single 
monomer by the usual 12-6 LJ expression [13]. Hence, 
the parameter e s weighs the monomer-surface (E sur f) and 
monomer-monomer (£?buik) interaction. Center-of-mass 
translation is restricted by the attractive substrate at 
z = and a sufficiently distant steric wall at z = L z . In 
our microcanonical analysis it will become clear how L z 
influences the results, however, the effect on the canonical 
data is small if L z exceeds the extension of the polymer. 
As long as not mentioned otherwise, the ratio N/L z is 
kept constant (L z = 37V). We always employ natural 
units (&b = 1). 

Energetic and Structural Quantities. To describe the 
canonical equilibrium behavior, we use the canonical ex- 
pectation values and thermal fluctuations of the following 
quantities: energy and specific heat, cy, the radius of gy- 
ration, (iigyr), as a measure for the extension of the poly- 
mer, and its tensor components parallel and perpendicu- 
lar to the surface, (R\\) and (R±), with i?g yr = -Rjj 
The components are of interest due to the structural 
anisotropy introduced by the substrate. Other useful 
quantities are the distance of the center-of-mass of the 
polymer to the surface, (z cm ), and the mean number of 
monomers docked to the surface. 

Multicanonical Sampling. The density of states g(E) 
encodes all information regarding the phase behavior of 
the system rendering its precise estimation extremely 
helpful. This requires the application of sophisticated 
Monte Carlo methods. In this work, we have performed 
multicanonical simulations [lOj. The idea is to increase 
the sampling rate of conformations being little favored in 
the free-energy landscape by performing a random walk 
in energy space. This is achieved by introducing suitable 
multicanonical weights W muca (E) ~ g~ 1 (E) to sample 
conformations X according to a transition probability 

w(X -> X') = mm[W muca (E(K'))/W muca (E(K)), 1]. 

(2) 

As the weights W muca (E), i.e. g(E), are unknown a pri- 
ori, they are determined iteratively until the energy his- 
togram is constant up to a variation of about 10% in the 
desired energy range. An efficient, error- weighted multi- 
canonical recursion is described in Ref. [14f . 



III. PSEUDO-PHASE DIAGRAM 

To construct the conformational pseudo-phase dia- 
gram, multicanonical simulations [lOj for 51 different sur- 
face attraction strengths e s £ [0, 5] were performed for a 
chain with N = 20 monomers. These data can now be 
reweighted to arbitrary temperature, but since it turns 
out that the interval T £ (0,3] is the most interesting 
one, we restrict ourselves to this range here. Each simu- 
lation consisted of 10 8 sweeps and was performed with at 
least two different initializations. The final pseudo-phase 
diagram is shown in Fig. [T] and representative conforma- 
tions are given in Fig. [2j The blue bands indicate the ap- 



2 




1 2 3 4 5 



FIG. 1: Pseudo-phase diagram of the 20mer (for details, see 
text). Different representative conformations are shown in 
Fig. [2] Along the lines of constant e s a microcanonical anal- 
ysis has been performed, see section HVl 

proximate phase boundaries that have some uncertainty 
because the peaks of the fluctuation of canonical expec- 
tation values do not coincide for finite systems. It should 
be stressed that due to the finite chain lengths all phases 
and transitions here are not phases in the strict ther- 
modynamic sense. Nevertheless, a reasonable picture of 
polymer adsorption behavior is obtained and most of the 
phases are believed to still exist for longer chains. Here, 
only some representative observables used for the con- 
struction will be discussed. For more details, see Ref. |{|. 

Energetic Fluctuations. Although the energy varies 
smoothly with T and e s , two transitions can be iden- 
tified as ridges in the profile of the specific heat: The 
adsorption transition separating desorbed and adsorbed 
conformations and a freezing transition at low tempera- 
tures. Near T = 0.25, cy exhibits a pronounced peak in- 
dependently of e s . The crystalline shape of the structures 
below this peak additionally confirms its nature as freez- 
ing transition. However, to identify different crystalline 
shapes, a closer look at the conformational quantities is 
needed. 

Structural Fluctuations. The average radius of gy- 
ration (Rgyr) reveals that the most compact conforma- 
tions dominate at low T and low e s . It establishes the 
phase boundaries between DE (desorbed expanded) and 
DG (desorbed globular) and between AE2 (adsorbed ex- 
tended; not flat on substrate) and AG (adsorbed globu- 
lar) and confirms the freezing transition, but the adsorp- 
tion transition is not prominently signaled by (Rgyr)- Its 
tensor components give additional information. For ex- 
ample, for e s > 3.4, (R±) vanishes at low T, whereas 
(R\\) attains low values at lower e s . Small values of (R±) 
correspond to conformations spread out flat on the sur- 
face, with associated pseudo-phases AC1 (adsorbed com- 
pact; flat) and AE1 (adsorbed expanded; flat), separated 
by the freezing transition. The most pronounced transi- 
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FIG. 2: Representative examples of conformations for the 
20mer in the different regions of the T-e s pseudo-phase dia- 
gram in Fig. [T] DE, DG, and DC represent desorbed "phases" . 
In regions AE1, AE2, AC1, AG, AC2a, and AC2b, conforma- 
tions are favorably adsorbed. 




FIG. 3: d (R±) /dT of the 20mer. 



tion is the strong layering transition at e s ^3.4 that sep- 
arates regions of planar conformations (AC1, AE1) from 
the region of double-layer structures (AC2b) and ad- 
sorbed globules (AG), below and above the freezing tran- 
sition, respectively. This sharp energetical transition can, 
e.g., be nicely seen in d (R±) /dT in Fig. [3l Although for 
the considered short chains no higher-layer structures are 
observed, (R\\,±) indicate some activity for lower e s . For 



N = 20, e s w 1.4 is the lowest attraction strength, where 
stable double-layer conformations are found. What fol- 
lows is a low-temperature subphase of surface attached 
compact conformations (AC2a). These structures occur 
if the surface attraction is not strong enough to induce 
the formation of compact layers. Structures here are sub- 
ject to quite strong finite-size effects. Raising the temper- 
ature above the freezing transition starting in the AC2 
regions, polymers adopt the adsorbed, globular, but un- 
structured conformations of the AG phase. This pseudo- 
phase has been first conjectured from short exact enu- 
meration studies of 2D polymers in poor solvent [15IL but 
was also found in lattice-polymer simulation studies[6|, H[ • 
At even higher T, two scenarios can be distinguished de- 
pending on the relative strengths of £t>uik and E sur f • For 
low e s , the polymer first desorbs (from AG to DG) and 
expands at even higher temperatures (from DG to DE). 
For larger e s , the polymer expands while it is still ad- 
sorbed (from AG to AE2) and desorbs at higher T (from 
AE2 to DE). The remaining observables confirm the pic- 
ture sketched so far. The center-of-mass distance to the 
surface (z cm ) gives a clear signal of the adsorption transi- 
tion, whose location is well described by T a d s oc e s . Since 
at higher T the stronger thermal fluctuations are more 
likely to overcome the surface attraction, this is intuitive. 
The mean number of surface contacts supports the ob- 
served layering. 

It is clear that in particular in the compact pseudo- 
phases the structural behavior of the studied small chains 
is affected by finite-size effects. However, especially 
at high temperatures, the pseudo-phase diagram con- 
structed here corresponds quite well with a similar lattice 
study Q with the advantage of not suffering from lattice 
artifacts. 



IV. THE ADSORPTION TRANSITION 
REVISITED MICROCANONICALLY 

We now concentrate on the adsorption transition and 
look at it from another perspective: the microcanonical 
one. This approach has already proven quite useful for 
first-order-like structural transitions such as molecular 
aggregation processes [HI 03 an d protein folding [l8l[l9j. 
For more details on this work see Ref. [20| . 

The central quantity is the density of states g(E) or the 
microcanonical entropy denned as S (E) = In g(E). Here, 
we normalize everything by the number of monomers and 
use 

s(e) = N- 1 lng(e), (3) 

with e = E/N . In contrast to canonical (NVT) statis- 
tics, where T is an externally fixed control parameter, in 
the microcanonical (NVE) ensemble it is derived from 
the entropy, T(e) = [ds(e) / de\^ v . There are cases for 
finite systems, where s(e) is a convex function in a tran- 
sition regime. A consequence is that with increasing sys- 
tem energy the temperature decreases. This is true as 
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FIG. 4: Microcanonical entropy s(e) (up to a constant) for 
a 20mer at e s — 5, the Gibbs hull 7i s (e), and the difference 
As(e) = 7is(e) — s(e) versus energy per monomer e. The local 
maximum of As(e), called surface entropy As sur f, defines the 
energy of phase separation. The latent heat Aq is defined as 
the energy being necessary to cross the transition region at 
the transition temperature T a a s - 

long as the surface-to- volume ratio is large enough to sup- 
press a concave increase of s(e). In such a case, the ener- 
getic separation of the two distinct phases is sufficiently 
large to establish a kinetic barrier. This regards all first- 
order phase transitions and two-state systems, but also 
transitions, where phase coexistence is completely absent 
in the thermodynamic limit, but not for the finite system. 
The latter is the case here: The adsorption transition of 
flexible polymers to an attractive substrate is known to 
be continuous in the thermodynamic limit. However, as 
we will show here, the adsorption of finite non-grafted 
polymers exhibits signals of a first-order transition which 
vanish in the thermodynamic limit. 

Exemplified for a 20mer and e s = 5, we have plotted 
in Fig. [4] the microcanonical entropy s(e). It shows the 
characteristic features of a transition with phase coexis- 
tence in a small system. For energies right below e a d s , 
the system is in the adsorbed phase AE2 (cf. Fig. [1]), for 
e a ds < e < edes, the system is in the transition region, 
where s(e) is convex. One can construct the Gibbs hull 

n s (e) = s(e ads ) + e(ds/de) e = eads (4) 

as the tangent that touches s(e a d s ) and s(ea e s), whose 
inverse slope T a d s = (dH s / de) -1 is the microcanoni- 
cal definition of the adsorption temperature. However, 
the transition rather spans a region of temperatures like 
the fluctuation maxima do in the canonical ensemble. 
Hence, this adsorption temperature definition is not the 
only one possible. A unique transition point only ex- 
ists in the thermodynamic limit. Nevertheless, not only 
for systems, where the thermodynamic limit is unreach- 
able [2l| in principle such as for proteins, it is worthwhile 
to understand the behavior of such a quantity. For a fur- 
ther analysis, we also use the surface (or inter facial) en- 
tropy, representing the entropic barrier of the transition, 




FIG. 5: s(e) and its fraction for adsorbed conformations 
Sads(e) at various surface attraction strengths e s = 0, 1, . . . , 6 
for a 20mer [for all e a , the fraction for desorbed structures, 
Sdes(e) resembles s(e) for e s = 0]. 

As surf = max{As(e) = H s (e) - s(e) | e ads < e < e des } 
and the latent heat, Aq = e^es — e a d s - Before we show 
for the adsorption transition that Aq decreases with TV, 
we first investigate the origin of the phase separation for 
finite chains. 

Dependence on the Surface Attraction Strength. In 
Fig. [5j s(e) is shown for a 20mer and different e s , where 
here also e s = 6 has been included. Since the high-energy 
regime is dominated by desorbed conformations, s(e) is 
hardly affected by e s here, while the low-energy tail in- 
creases significantly with e s . Thus, one can split the 
density of states into contributions of desorbed and ad- 
sorbed conformations, #des(e) and g a ds(e), respectively, 
such that g(e) = gdes(e) + #ads(e) and s des ,ads(e) = 
TV -1 In 

#des,ads( e )- We consider the polymer to be ad- 
sorbed if E sur f < — 0.1e s iV. This choice includes all 
polymers that are reasonably strongly adsorbed in terms 
of energy and works for all e s . Moreover, qualitative 
features do not depend sensitively on the choice and this 
devision is only employed to get a qualitative picture, not 
to extract Aq or As sur f. Since both, s a d s (e) and Sdes(e), 
are concave in the whole energy range of the transition, 
the convex entropic monotony can only occur in the re- 
gion, where adsorbed and desorbed conformations have 
similar entropic weight. 

Performing the Gibbs construction as in Fig. [4] to ex- 
tract As(e), one sees [20] that for e s ^ 2 the transition 
appears to be first-order like (Aq = e^es — e a d s > 0) 
for a finite, non-grafted chain. For e s ^ 2, the Gibbs 
construction is no longer meaningful in absence of a con- 
vex regime in s(e), indicating a second-order phase tran- 
sitions (Aq = 0). Referring to the phase diagram in 
Fig. [TJ the adsorption transition thus seems to become 
first-order-like at that point, where it falls together with 
the O-transition (e s « 1.8, T « 1.3). This is also signaled 
by the saddle point of the corresponding T~ 1 (e) curve. 
For larger e s , phase coexistence gets apparent between 
DE and AE2. Here, As sur f and Aq increase with e s and 
trivially diverge for e s — >• oo. Also the first-order-like 
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features of T _1 (e) increase and the adsorption tempera- 
tures T a ds depend roughly linearly on e s , as was already 
suggested by the canonical data. 

Chain-length Dependence. Since the adsorption tran- 
sition is expected to be of second order in the thermo- 
dynamic limit [5|, first-order signatures found for the fi- 
nite system between DE and AE2 must disappear for 
N —> oo. Indeed, our data for N up to N = 150 support 
a power-law scaling of the latent heat, Aq ~ N~ k q, with 
K q ~ 0.35 — 0.40, which clearly suggests liniTv^oo Aq = 0, 
confirming this expectation. 

Variation of the Box Size. After noticing that there 
is a considerable influence of the simulation box size on 
s(e), we also investigated this effect. To this end, simu- 
lations with fixed e s = 5 and chain length N = 20 were 
performed for different L z = 20, 30, . . . , 150. Because 
the number of adsorbed conformations cannot depend 
on L z , the unknown additive constants to s(e), s a d s (e), 
and Sdes(e) were chosen such that s a d s (e) coincides for 
all L z . With this choice, Sdes(e) increases with the loga- 
rithm of L z , like it should be the case for the translational 
entropy. Consequently, both, the surface entropy As sur f 
and the latent heat Aq increase with L z . Note, that in 
the case of a grafted polymer, effectively corresponding 
to a small L Z1 we did not observe any convex intruder in 
the microcanonical entropy. 

V. SUMMARY 

In this work, we have used two approaches to describe 
the behavior of a single homopolymer near an attractive 
substrate. 

First, in analyses of canonical expectation values of 
several energetic and structural quantities and their ther- 
mal fluctuations for a chain with 20 monomers, confor- 
mational phases and phase boundaries in the pseudo- 



phase diagram versus temperature and surface attrac- 
tion strength were identified. Our chosen simulational 
method was the multicanonical Monte Carlo technique. 
Although the computational expense to accurately ex- 
plore such a broad parameter range restricted us to rather 
short chains, for the majority of pseudo-phases, in par- 
ticular those that are assumed to be relevant in the ther- 
modynamic limit, we find a nice qualitative coincidence 
with similar lattice studies. Then, we complemented the 
picture by focusing on the adsorption transition micro- 
canonically. For short polymers, the microcanonical en- 
tropy revealed that at the adsorption transition adsorbed 
and desorbed conformations coexist, corresponding to a 
first-order character of this transition for short polymers. 
We have studied how the character of this transition de- 
pends on surface attraction strength, chain length, and 
concentration. 

Altogether, our study has shown the usefulness of a 
combined approach of the microcanonical and canonical 
ensemble in understanding the conformational behavior 
of finite systems. 
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